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Abstract 

We study long-range correlation functions of the rectangular Ising lattice with cyclic 
boundary conditions. Specifically, we consider the situation in which two spins are on the 
same column, and at least one spin is on or near free boundaries. The low-temperature 
series expansions of the correlation functions are presented when the spin-spin couplings 
are the same in both directions. The exact correlation functions can be obtained by D log 
Pade for the cases with simple algebraic resultant expressions. The present results show 
that as the two spins are infinitely far from each other, the correlation function is equal to 
the product of the row magnetizations of the corresponding spins as expected. In terms 
of low-temperature series expansions, the approach of this m-th row correlation function 
to the bulk correlation function for increasing m can be understood from the observation 
that the dominant terms of their series expansions are the same successively in the above 
two correlation functions. The number of these dominant terms increases monotonically 
as m increases. 
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1 Introduction 



The Ising model is one of the simplest models in the study of critical phenomena, and has been 
investigated by many authors since Onsager's celebrated work on the free energy of the rectan- 
gular lattice without magnetic field [1] , followed by Yang's derivation of the spontaneous mag- 
netization [2] which gave the critical exponent (3 — 1/8. The correlation function < cro,oO"m,n > 
of the rectangular lattice has also been well studied, and can be expressed in terms of Toeplitz 
dctcrminantal form with dimension 2m + 2n [3, 4, 5] and a "generalized Wronskian" form with 
dimension m [6]. The most efficient way to calculate the bulk correlation functions is given in 
Refs. [7, 8] and the references therein. McCoy and Wu considered the rectangular lattice with 
cyclic boundary conditions, and only the spins on a boundary row interact with a magnetic 
field [9]. In this case, the critical exponent of the boundary spontaneous magnetization was 
found to be 1/2, and the magnetization in the m-th row was expressed by a determinant with 
dimension m. However, the approach of the boundary spontaneous magnetization to the bulk 
result was unsolved. As the square of the spontaneous magnetization is equal to the long-range 
correlation function, this problem is equivalent to the approach of the boundary correlation 
function to the bulk correlation function. It is the purpose of this paper to study this approach 
by the long-range correlation functions with at least one spin on or near the free boundaries. 
In section 2, we give the known and expected results of the correlation functions together with 
the basic definition. The calculation and expressions of the correlation functions between two 
spins on a column will be sketched in section 3. Since the correlation functions do not have a 
simple closed form in general, we content ourselves with the expression of the correlation func- 
tions in terms of low-temperature series expansions, and indicate the approach of the boundary 
correlation function to the bulk correlation function in section 4. Conclusions will be given in 
Section 5. 



2 Rectangular Ising lattice and the correlation functions 



We mainly consider here a rectangular lattice with cyclic boundary conditions imposed in the 
horizontal direction, and free boundary conditions in the vertical direction without magnetic 
field. The numbers of rows and columns are denoted by M. and jV, respectively. The Hamilto- 
nian of the system is 



where Ji and J2 arc the horizontal and vertical spin-spin interactions between nearest-neighbor 
spins, and o-^^^+i = (jj^i. Each a assumes the value +1 or —1. 

In the two-dimensional limit, i.e., M. ^ 00 and A/" — > 00, the solution of the spontaneous 
magnetization is well-known [2, 3, 10], namely the bulk correlation function < am,ic^m,n > 
between two spins on the same row in the limits n— >oo,j\/'— n— >oo and m— >oo,A1— m— >oo 
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where i^i = pJi,K2 = (3J2 and Zi = tanhKi, z^ = tanhi^2- 

The boundary correlation function < (J\^\a\^n > was considered in [9]. 
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< cri,i^i,n >= 



<i 



Azizl 

cosh(2X2) - coth(2Ki) 



cosh (2X2) - 1 

The second-row correlation function < cr2,icr2,n > in the limit n — > 00, — n — > 00 is [9] 
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where a is the larger root in magnitude of the equation 

{l + zl){l + zl)-z^{l-zl){e 



)- Z2{\- z{){a-Va-') =0 
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Let us consider the long-range correlation function between two spins on the same col- 
umn. In the limit that the distance between these spins is infinite, the correlation function 
should be the product of the magnetization of each spin. Therefore, the bulk correlation func- 
tion < o"m',iCm,i > with m' < m in the limits m — m' — > cxo, m' — > 00, — m — > 00 
should be < o-m',iO'm,i >= ^s^buik- Similarly, we expect the following results: the boundary- 
boundary correlation function Xmij^^ao < o'i,iO'a^,i >= the second row-second row cor- 
relation function lim^vi^oo < c"2,iC"a^-i,i >= M^,2j ^^d the boundary-second row correlation 
function limji^_^^ < a2,icrM,i >— -^s,i-^s,2, etc. If one of the spins is far from the free 
boundaries, say on the middle row, then we expect that the boundary-middle correlation func- 
tion lim^K^oo < (^i,i<^[M/2],i >= ^s,i^s,buik and the second row-middle correlation function 
lim^^oo < cr2,icr[M/2],i >= Ms,2Ms,buik, etc., where M^^buik, M^^i, Ms,2 are given in eqns. (2.2) to 
(2.4) and [g] denotes the integral part of g. Whether A4 is even or odd is not important here 
since we take the limit Al — > 00. 



3 Calculation for the correlation function between two 
spins on the same column 

Wc calculate correlation functions between two spins on the same column by the standard 
Pfaffian approach. Since this method is available in the literature [3, 4, 9, 11] for decades. 
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the derivation of our results should be brief. As is well-known, the rectangular Ising lattice is 
equivalent to the dimer problem on the "bathroom tile" lattice. The square of the partition 
function is given by 

= (2coshXi)^-^(coshK2)^(^-^) det^ (3.1) 

where A is an antisymmetrical matrix, and some of its elements in the limit — > oo will be 
given explicitly in eqns. (3.5) to (3.8). The square of the correlation function between cr^/^i 
and (7m,\ with m' < m can be shown to be 



< <Tm',i<Tm,i zf"'^ det(r^ + A'^) det 5 
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where in the standard notation 
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so that det (5 = (2:2 "'^—2;2)^*-™ '"'^and^ Ms the submatrix of ^ in the subspace {(m', 1);7, (m, l)z)}. 
Explicitly, we have 



< Crm',lCFm,l > = \} - ^2) ^ 
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(m' + l;m)DD 
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where wc use the abbreviations Z2 



Z2 



1-1 



and {j;j')xx' = A \j,l;j',l)xx' with 



X,X' G {U,D}. By a procedure similar to that in [9], the elements of the inverse matrix A ^ 
can be calculated. In the limit A/" — > 00, they are 

A~\j, k;j', k')DD = -A~^{3', k';j, k)DD 
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A-\j,k;f,k')uu = -A-\f,k';j,k)uu 

A-\j, k;j', k')uD = -A-^ij', k';j, k)Du 

1 r dee^e(k-k'^''''~'^^ - - 11 + ^le^Tf ) [1 + a-'(M^^ni + «-2(-M-.)(^)2] 



(3.7) 

A-\j, k;j', k')Du = -A-\3\ k';j, k)uD 

27r7_/^" (a-i - a)z2(l - zf) i + «-2Ai(^)2 ^-^-^^ 

where 1 < j' < j for the first three functions, and 1 < j' < j for the last one. c and c' are 
defined as 



2 _L ^2 _ ^2 \ , / 1 z| + _ 52 



where a = 2i;2isin6'|l + zie'^l'^, 6 = (1 - zl)\l + zie'^\-'^, and A = |1 + zie'^\'''^Z2{l - zj)a, 
X'^\l + Zie''\-^Z2{l-zl)a-\ 

Because A~^{j, k;j', k)ciD = A^^{j, k;j', k)uu = and the sign of < (Jm',iO'm,i > should be 
chosen so that < (Jm' ,i'^m,i 1 as T — > 0, we obtain the following general expression of the 
correlation function between two spins on the same column, 

Z2 — Z2A~^{m';m' + l)uD —Z2A~^{m';m' + 2)uD ■■■ —Z2A~^{m';m)uD 
-Z2A'^{m' + l;m' + l)uD Z2 - Z2A~^{m' + l;m' + 2)ud ■■■ -Z2A~^{m' + l;m)uD 

-Z2A~^{m-l;m' + l)uD -Z2A~^{m - l]m' -\-2)ud ■■■ Z2 - Z2A~^{m - l;m)uD 

(3.10) 

where we define 

-Z2A-\j;j')uD = -(1 - zl)A-\j, 1;/, 1)ud 
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for 1 < j' < j 



-L ~r Z\ 6 I 



(3.11) 
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4 Correlation functions in terms of low- temperature se- 
ries expansions 

The correlation function < cTm',iC"m,i > given in eqn. (3.10) is expressed by the determinant 
of a (m — m') x (m — m') matrix. As we are interested in long-range correlation functions, 
the size of the matrix approaches infinity. Instead of performing the integrations of all the 
matrix elements directly, we expand the integrands in terms of low-temperature expansions. 
For simplicity, we consider Ji = J2 = J ■, and use the low-temperature variable z = e~'^^'^ . It is 
expected that the first few terms of the series for a finite value oi m — m' become the same as 
in the infinite case and that number of correct terms would increase as the difference m — m' 
increases. If the exact result has such a simple expression as in eqns. (2.2) or (2.3), it can be 
obtained by D log Fade analysis [12] with the first certain terms known. We shall check this 
procedure on the bulk correlation function first. Then we apply it to other cases in this section. 



4.1 Bulk case 

For the bulk correlation function, we take the limits m — m' — > 00, m' — > 00, M. — m — > 

00. The last two limits will not be written out explicitly in this subsection. In these limits, 
the last three factors in eqns. (3.11) and (3.12) are equal to 1. Using the notations /j-j' 
and fj_j, for the corresponding matrix elements (1 — 2;|).4~^(j', 1; j, 1)ud for I < j' < j and 
(1 — Z2)A~^{j, 1;/, l)uD for 1 < j' < j, respectively, the bulk correlation function is given by 



lim < am',icrm,i > 



m—m—*oo 



lim 

m—m'—*oo 



Z2 — fl —f2 —fs • • • —fm-m' 

~fo ^2 — fl — /2 • • • —fm-m'-l 

~fi ~/o ^2-/1 ■ ■ : 

■■ '■■ ■■■ ■■•^ -/2 

~fm-m'-2 ~fm-m'-3 ' ' ' ~fo ^2 — fl 

(4.1.1) 

This is the traditional Toeplitz determinant, and the matrix corresponds to the infinite-size 
center part of the matrix for limx_^oo < cri,iCM,i > ^ shown in the next subsection. We 
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consider here the £x£ submatrix of the above matrix by keeping the first £ columns and the first ^ 
rows where I — m—m' is finite. Then, we define the correlation function lim_A4^oo < crm',iO'm,i >e 
by the determinant of that Ixl submatrix. By expanding the integrand of each matrix element 
in the power series of 2;, by performing the integrations and by computing the determinant, we 
obtain 



lim 
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1-4^^ 


- 12/ 
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>3 
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6Az^ - 


268z^° 


- 1152z^^ 
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< CFm' 


,lCm,l 
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16z^- 
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272z^° 


- 1224z^^ 


- 5680^^^ + ... 


lim 
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,icrm,i 
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= 1 - 
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16/- 


64;^^ - 


272^1° 


- 1228;^^^ 


- 5788/^ + ... 


lim 

A^— >oo 


< CTm' 


,l'^m,l 


>6 


= 1 - 


4^^- 


16/- 


64^^ - 


272^^° 


- 1228^^2 


- 5792z^^ + ... (4.1.2) 



Compared with the exact bulk correlation function 

1/4 

(l-/)4_ 

= 1 - 4/ - 16/ - 64/ - 272z^° - 1228^^^ - 5792z^^ 



liin < am'lCTmA > = 



m—m' ^00 



^ 16/ 



-28192^^^ - 140448^^^ - 712276/° + ... (4.1.3) 

We find that the series of \\m.j^^oo < 0"m',iCm,i >f is correct up to the z'^^^'^ term. Because the 
coordination number of non-boundary spins is 4, the power of z in the series expansion terms 
for the bulk correlation function is always even. In this case, u = z'^ is used as the expansion 
variable traditionally. Notice that this is no long true for the other cases we consider in the 
rest of this section. Since 

-^ln( lim < (T^',ia^,i >) = 2\n^^a ^ 2^ ' (^-l-^) 

du Vm-m' -»oo / [1 — U^)[l — OU + U^) 

one only needs terms up to 2;^^ to obtain the exact result by D log pade because only 5 unknown 
parameters in this equation need to be determined. 



4.2 Boundary-boundary case 

In order to obtain the boundary-boundary correlation function, we substitute m' = 1 and 
m = Ai in cqn. (3.10). Wc show here explicitly the power series of the correlation function 
< cn.io'x,! >M which is given by the determinant of a {M. — 1) x {A4 — 1) matrix (3.10): 

< (Ji,i(J2,i >2 = 1 -4/ - 8/ -4/ + 48/ + 184/ + 272/ + ... 
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< (^1,10-3,1 >3 


= 1 
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- 16/ 


-24/ 


- 16^^ 


+ 152^' 


' + ... 


< ('■1,1(74,1 >4 


= 1 


/( 3 

— Az - 


o 4 
- 82; 


— loz 


r\ A 6 

— 242; 


— 362; 
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— 482; 


+ ... 


< (^1,1(^5,1 >5 


= 1 


-4/ - 


-8/ 


- 16/ 


-24/ 


-36/ 


-48/ 


-64^' 


< (^1,1(^6,1 >6 


= 1 


-4/ - 


-8/ 


- 16/ 


-24/ 


-36/ 


-48/ 


-64^' 








+ ... 













(4.2.1) 

Compared with the following exact result in the limit — > 00 

(l-2z -/)(! + /) 
hm < ai,iaM,i > = r^j^ — ^ 

M^oo (1 — 2;)'^(1 + Z) 

= l-Az^ - 8z^ - 16z^ - 24:z^ - 36z'^ - A8z^ - 6Az^ 

-80z^° - 100/^ - 120/^ - 144^^3 - 168/^ + ... (4.2.2) 

we observe that < (Ji,i(Tm,i >m is correct to the 2;^-^ term. Now one only needs terms up to 
z^^ to obtain the exact result by D log Fade because 

d ^ / ^. \ -8/(3 + 2^ + /^ 
— m hm < (Ti 1 cT/u 1 > = tz tttz z 



"M rr- < (^i,i(^M,i > ) = — ^ ^ • (4.2.3) 

dz \M-^oo / [1 — z^){l — 2z — Z^) 

4.3 Second row-second row and higher row cases 

For the second row-second row correlation function, < o'2,iO'x-i,i >, its matrix corresponds 
to the matrix for < cri,i(Tx,i > with the first and the last columns and rows removed. Again 
consider finite M. so that < 02,1(7 m-i,i >m is the determinant of the {M.—3) x {M. — 3) matrix. 
The cases = 4, 5 and 6 are given as 



< (^2,1(^3,1 >4 = 


1 


-4/- 


4^5 


-24/ 


- 52^^ - 


112/ + ... 




< Cr2,lC4,l >5 = 


1 


-4/- 


Az^ 


-28/ 


- 56^^ - 


156/ - 296/ 


- 576^^° + 


< Cr2,l05,l >6 = 


1 


-4/- 


4z^ 


-28/ 


- 56^^ - 


160/ - 300^^ 


- 648^^° 



-1124/^- 2076/^ + ... (4.3.1) 
These series expansions approach the expected result for — > 00 which is 

lim < (72,1(7^-1,1 > = 1 - 4/ - 4/ - 28/ - 56/ - 160/ - 300/ - 652z^° - 1128z^^ 

-2188^^^ - 3500^^^ - 6340^^^ + ... (4.3.2) 
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Notice that there is no term in contrast with the boundary-boundary case in eqn. (4.2.2), 
and the coefficient of z'^ term is the same as the coefficient of z^ term for the bulk case in eqn. 
(4.1.3). 

In general, the matrix for < crm+i,ic>i-m,i > corresponds to the matrix for < ai^i(TM,i > 
with the first and the last m columns and rows removed, and we expect that lim^^oo < 
Cm+i,iCAi-m,i > is the same as the square of the (m + l)-th row magnetization. However, there 
may be no simple expressions for these correlation functions, and these can not be obtained by 
D log Pade analysis. 

As the next example, the matrix for the third row-third row correlation function is the 
submatrix of the matrix for < (Tia'^m,! > with the first and the last two columns and rows 
removed. The first two series expansions of < 0-3 i(Tyv(-2,i >m with finite Ai are 

< (73,iO-4,i >6 = 1 - 4z^ - 12/ - 4^^ - 56/ - 124/ - 432z^° + ... 

< ^^3,1^5,1 >7 = 1 - 4/ - 16/ - 4/ - 80/ - 128/ - 564^^° + ... (4.3.3) 

In this case, there is no z"^ and z^ terms, and the coefficients of z^ and z^ terms are the same as 
the corresponding coefficients for the bulk case in eqn. (4.1.3). Therefore, we expect that the 
number of such leading terms in the series expansion of limTn^oo < o'm+i,iO'M-m,i > as agree 
with the corresponding terms for the bulk case increases monotonically as m increases. 

4.4 Boundary-second row and boundary- higher row cases 

For the boundary-second row correlation function, < a2^i(JM,i its matrix corresponds to the 
matrix for < ai^iaM,i > with the first column and row removed. For finite M., the correlation 
function < (T2,i(Tx,i >m is given by the determinant of the {Ai — 2) x {M. — 2) matrix. The 
power series expansions for small values of A4 are given in the following: 



< cr2,icr3,i >3 


= 1 


-2/ 


-6/ 


-8/- 


- 18/ - 


-20/ + 56/ + ... 






< 0'2,lC^4,l >4 


= 1 


-2/ 


-6/ 


- IQz^ 


-28/ 


- 48^^ - 104/ + 






< (^2,10-5,1 >5 


= 1 


-2/ 


-6/ 


- 10/ 


-28/ 


- 50/ - 118/ - 


192/ 


- 362z^° + 


< CF2,l(y&,l >6 


= 1 


-2/ 


-6/ 


- 10/ 


-28/ 


- 50/ - 118/ - 


194/ 


- 380^^° 



-540/^- 928/2 + ... (4.4.1) 
These series expansions approach the expected result for ^ 00 which is 
lim < (72,i(7jw,i > = 1 - 2/ - 6/ - 10/ - 28/ - 50^^ - 118/ - 194/ - 380^^° - 542^^^ 

A^— >oo ' ' 

-950/2 _ 1142^13 _ 1960^14 + (4 4 2) 

Notice that the coefficient of z^ term is half the value of the coefficient of z^ term of the 
boundary-boundary case. 
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In general, the matrix for < (Jm'+i,i<^M-m,i > corresponds to the matrix for < ai^iaM,i > 
with the first m' columns and rows and the last m columns and rows removed, and we expect 
that lim^^oo < crm'+i,icrM-m,i > is the same as the product of the (m'+l)-th row magnetization 
and the (m + l)-th row magnetization. 

As another example, the matrix for the boundary-third row correlation function is the 
submatrix of the matrix for < cri^ia^vj,! > with the first two columns and rows removed. The 
first three series expansions of < lO";^ i >m with finite A4 are given by 

< 0-3,1(74,1 >4 = 1 - _ g^4 _ g^5 _ ^2z^ _ 4^7 _ ^q^8 ^ 

< (^3,1(^5,1 >5 = l-2z^ -6z^ -8z^ -22z^ -22z'^ -70z^ -SAz^ -260z^^ + ... 

< (^3,1(^6,1 >6 = l-2z^- 6z^ - 8z^ - 22z^ - 24/ - 84/ - 122/ - 420^^° - 712z^^ 

-m4z^^ + ... (4.4.3) 

In this case, the coefficients up to z^ terms are the same as the corresponding coefficients for 

Ms.i^s.buik given in eqn. (4.5.2). We expect that the number of such leading terms in the series 
expansion of limyn^oo < o'm+i.iO'M,! > as agree with the corresponding terms for Ms,iMs,buik 
increases monotonically as m increases. 



4.5 Boundary-middle case 

For the correlation function between one spin on a free boundary and the other spin on the 
middle row, substitute m' = 1 and m = [M./2] in eqn. (3.10). The matrix for < cri,i(T[x/2],i > 
corresponds to the upper-left quarter of the matrix for < cri,i(Tx,i >. In the limit M. — > oo, 
the two factors of eqns. (3.11) and (3.12) containing a to negative infinite powers are equal 
to 1. As before, lim^K^^oo < <^i.i(^[m/2],i >e denotes the determinant of the £ x £ submatrix of 
the matrix for lim_A4^oo < ^i,iO'[m/2],i > with first i columns and first £ rows. The expansion 
formulas for small £ are given as 

lim < (7i,i(7r^/2] 1 >i = 1 - 2/ - 6/ - 6/ - 12/ - 2/ - 6/ + 34/ + 36z^° + 162z^^ 

+170/2 + ... 

lim < (7i,i(JrA^/2],i >2 = 1 - 2/ - 6/ - 8/ - 22/ - 20^^ - 60/ - 20/ - 94^^° + 168z^^ 

+114/2 + ... 



^im < (Ti,iCT[^/2],i >3 = 1 - 2/ - 6/ - 8/ - 22/ - 22^^ - 74/ - 58/ - 254^^° - UOz 

-728^^2 ^ 

lim < (Ti,iCT[A^/2],i >4 = 1 - 2/ - 6/ - 8/ - 22/ - 22^^ - 74/ - 60/ - 272^^° - 176^ 
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-1090^^2 + ... 

lim < (7i i(7rjw/2] 1 >5 = l-2z^ - 6^^ - 8^^ - 22z^ - 22z'^ - 7Az^ - 60z^ - 272^^° - 178^^^ 

jM— »oo ' 

-1112^^2 + ... (4.5.1) 
which should be compared with the exact result, 

1/8 



lim < (7i,i(7[>i/2],i > = 



;i-2z-^2)(l + ^2) 
3? 



16z^ 



2\A 



= l-2z^- 6/ - 8^^ - 22/ - 22^^ - 74/ - 60/ - 272^^° - 178^^^ 
-1112^^2 - 580/^ - 4918/^ - 2018^^^ + ... (4.5.2) 

By D log Pade and with Zc — \f2 — 1, one only needs terms up to z^'^ to obtain the exact result 
because 

d,(y \ -2/(3 + 12^ + 2/-/) 

Tz ^Hi^^c^ < "^'^"[^/^i'^ n = (i-/)(i-2.-/)(i+2.-/) ■ ^'-'-'^ 

4.6 Second row- middle and higher row-middle cases 

For the second row-middle correlation function, < (T2,\(T[m/2],i >, its matrix corresponds to 
the matrix for < (Ji^ia[M/2],i > with the first column and row removed. Now we calculate the 
quantity lim^^oo < c'"2,iO'[yVi/2],i >e which is defined by the determinant of the i x i submatrix 
of the matrix for lim^K^oo < C2,iO'[a^/2],i > with the first i columns and the first i rows, we 
calculate 

J^im < (T2,i(T[A^/2],i >i = 1 - 4/ - 2/ - 18/ - 26^^ - 74/ - 94/ - 166^^° - 86^^^ 

-40/^ + ... 

lim < (72,i(7rjw/2],i >2 = 1 - 4/ - 2/ - 22/ - 28/ - 108/ - 148/ - 402z^° - 484^^^ 

-1082^^2 ^ 

lim < (72,i(7rjw/2],i >3 = 1 - 4/ - 2/ - 22/ - 28/ - 112/ - 150/ - 460z^° - 574z^^ 

-1688^^2 + ... 

lim < (72,iq>t/2],i >4 = 1 - 4/ - 2/ - 22/ - 28/ - 112/ - 150/ - 464z^° - 576^^^ 

-1778/^ + ... (4.6.1) 
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These series expansions approach the expected result which is 

Jim < (72,i(7[jK/2],i > = 1 - 4^^ - - 22z^ - 28^^ - 112^^ - IhQz^ - 464z^° - 576^^^ 

-1782^^2 - 2QUz^^ - 7054;^^^ - 7240^^^ + ... (4.6.2) 

Notice that there is no z'^ term in contract with the boundary-middle case in eqn. (4.5.2), and 
the coefficient of z"^ term is the same as the coefficient of z^ term for the bulk case in eqn. 
(4.1.3). 

In general, the matrix for < (Jm+i,iO^[M/2\,i > corresponds to the matrix for < o'i_i(T[_/vi/2],i > 
with the first m columns and rows removed, and we expect that lim^_>oo < crm+i,iC[M/2],i > is 
the same as the product of the bulk magnetization and the (m + l)-th row magnetization. 

As the next example, the matrix for the third row-middle correlation function is the subma- 
trix of the matrix for \imM~*oo < <7i,i0"[x/2],i > with the first two columns and rows removed. 
The first three series expansions of lim_yvi^oo < o'3,iO"[x/2],i >e are given by 

lim < (73,i(7[^/2],i >i = 1 - 4^^ - 12z^ - 2z^ - 46^^ - 62z^ - 276z^° - 518z^^ 

>oo ' ' 

-1278;2^2 + 



< (T3,i(7[^/2],i >2 = 1 - 4^^ - 16z^ - 2z'^ - 70z^ - 64:Z^ - 394^^° - 624z^^ 

-2052^^2 ^ 

lim < as iaiM/2],i >3 = 1 - 4/ - 16z^ - 2/ - 7Az^ - Uz^ - 438z^° - 626z^^ 

A^— >oo ' 

-2362z^^ + ... (4.6.3) 

In this case, no z^ and z^ terms appear, and the coefficients of z'^ and z^ terms are the same as 
the corresponding coefficients for the bulk case in eqn. (4.1.3). Therefore, we expect that the 
number of such leading terms in the scries expansion of limj^_>oo < crm+i,iC[A^/2],i > as agree 
with the corresponding terms for the bulk case increases again monotonically as m increases. 



5 Conclusions 

In conclusion, long-range correlation functions between two spins on the same column of the 
rectangular Ising model with cyclic boundary conditions were studied. The low-temperature se- 
ries expansions of the correlation functions were presented when the spin-spin couplings are the 
same in both directions. We confirmed that as the two spins are infinitely far from each other, 
the correlation function is equal to the product of the row magnetizations of the corresponding 
spins. The approach of this m-th row correlation function to the bulk correlation function for 
increasing m could be understood using low-temperature series expansions. Namely the dom- 
inant terms of their series expansions are the same and the number of these terms increases 
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monotonically as m increases. 
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